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SECTION  I 


INTRODUCTION 

One  failure  mode  experienced  with  gears  is  that  of  scoring.  Evidence  of 
this  failure  mode  is  the  existence  of  scratches  extending  outward  along  the 
working  surfaces  of  gear  teeth,  usually  from  the  pitch  line  to  the  tooth  tip 
(see  Fig  1).  (Occasionally  scoring  extends  below  the  pitch  line  to  the  tooth 
root.)  Continued  operation  when  the  surfaces  are  being  scored  will  lead  to 
significant  changes  to  the  tooth  working  surface  contour.  A  change  in  the 
surface  contour  may  be  expected,  at  the  very  least,  to  alter  noise  and 
vibration  levels. 

Implicit  with  the  change  in  surface  contour  is  the  generation  of  wear 
metal  particles.  These  particles  have  the  potential  to  further  damage  the 
operating  surfaces  of  the  scoring  gears,  other  gears  in  the  system,  bearings, 
seals,  etc.  The  final  result  could  be  catastrophic  failure. 


Destructive  Scoring:  Heavy  scoring  has  taken  place  above 
and  below  the  pitch  line,  leaving  the  material  at  the  pitch 
line.  As  a  result,  the  pitch  line  pits  away  as  it  attempts 
to  redistribute  the  load.  Usually  the  gear  cannot  correct 
itself  and  ultimately  fails. 


Fi gure  1 


Scuffed  Tooth  Surface 


(1) 
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It  is  generally  conceded  scoring  can  only  occur  when  the  lubricating  film 
is  no  longer  of  sufficient  thickness  to  prevent  metal-to-metal  contact  of 
mating  surfaces,  but  scoring  does  not  always  occur  in  the  absence  of  full  film 
lubrication.  When  conditions  are  conducive  to  the  formation  and  maintenance 
of  a  boundary  lubrication  film,  operation  can  continue.  In  some  cases  perfor¬ 
mance  may  even  be  improved  over  the  full  film  condition.  This  has  been  noted 
by  Godfrey,  Groszek,  and  others^'  .  Unfortunately,  the  phenomena  of 
boundary  lubrication  are  not  well  understood.  As  a  result,  designers  of 
lubricated  surfaces  tend  to  avoid  operation  in  the  boundary  regime,  yet  it 
cannot  be  avoided  entirely  for  during  starting  and  stopping  of  some  devices 
(such  as  gears,  cams,  bearings,  and  seals).  Operation  will  almost  always  occur 
in  the  boundary  lubricated  regime  unless  lubrication  is  by  hydrostatic  means. 
Of  course,  temporary  overloads  often  cause  operation  in  the  boundary  regime 
too. 
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SECTION  II 


GENERAL  DISCUSSION 
1 .  Definition  of  the  Objective 

The  field  of  boundary  lubrication  is  very  complex.  It  is  anticipated  that 
a  highly  interdisciplinary  approach  is  required  to  obtain  complete  under¬ 
standing  of  the  phenomena.  There  certainly  appears  to  be  considerable 
difficulty  in  understanding  the  details.  Understanding  seems  lacking  in 
areas  such  as:  film  formations  and  destruction;  surface  texture  effects; 
lubricant/metallurgy  interaction  effects;  and  operational  environment  effects. 
The  level  of  difficulty  in  understanding  boundary  lubrication  may  be  appreciated 
when  Beerbower's  depiction  of  "Regimes  and  Modes  of  Wear  on  Broken-In  Machine 
Parts"  is  viewed.  Beerbower  says  the  dashed  lines  in  this  figure 
(Ref:  Fig  2)  are  "both  arbitrary  and  ill  defined" In  this  figure: 


Where : 

V  = 

3 

Wear  volume,  cm 

d  = 

Sliding  distance,  cm 

Pm  = 

2 

Hardness,  kg/m 

W  = 

Load,  kg 

5  = 

Thickness  of  lubricant  film,  m 

D1  = 

Surface  roughness  of  body 

1 ,  ym 

°2  " 

Surface  roughness  of  body 

2,  ym 
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A  more  recent  report,  Lubricant/Metallurgy  Interaction  Effects  on  Turbine 
Engine  Lubricant  Load  Rating"^  emphasizes  gear  tooth  failure  by  scuffing; 
therefore,  a  considerable  portion  of  the  effort  focused  upon  boundary 
lubrication.  The  reported  effort  included  a  literature  review,  analysis,  and 
experiments.  The  primary  experimental  tool  used  was  the  Air  Force  Aero 
Propulsion  Laboratory's  Disc  Tester.  This  tool,  similar  to  others  in  the 
industry,  allows  simulation  of  the  sliding/rolling  contact  zone  of  a  lubri¬ 
cated  device  under  conditions  reasonably  well  controlled  and/or  amenable  to 
measurement. 

Typical  variables  of  interest,  controlled  or  measured  to  varying  degrees, 
i ncluded: 

a.  Sum  Velocity.  Literally  the  sum  of  surface  tangential  velocities  at 
the  line  of  contact. 

b.  Sliding  Velocity.  The  difference  of  the  surface's  tangential 
veloci ties . 

c.  Lubricant  Reactivity.  The  chemical  behavior  of  the  device's  surface 
(material)  with  the  lubricant. 

d.  Surface  Temperature.  The  average  of  the  mating  surfaces. 

e.  Conjunction  Temperature.  The  temperature  in  the  load  or  contact  zone 

f.  Surface  Topography.  The  roughness,  waviness,  and  lay  of  the  surface. 

g.  Lubricant  Supply  Temperature.  The  temperature  of  the  fluid  as  it  is 
discharged  from  the  nozzle. 

h.  Heat  transfer  characteristics  of  the  lubricant  and  the  device. 

It  was  not  completely  determined  how  the  above  listed  parameters  relate 
and  interact  to  allow  successful  operation  in  the  boundary  lubrication  regime 

It  is  also  interesting  to  note  that, after  consideration  of  the  state  of 
the  art  and  their  own  work,  the  authors  of  Lubricant/Metallurgy  Interaction 
Effects  on  Turbine  Engine  Load  Rating  concluded:  "The  basic  mechanisms  of 
scuffing,  pitting,  and  wear  under  steady  operating  conditions  such  as  in  the 
sliding/rolling  disc  system  are  far  from  understood." 

A  review  of  the  reports  previously  identified  and  other  literature  will 
confirm  this  conclusion.  The  conclusion  is  further  reinforced  by  the  various 
differences  of  opinion  on  the  meaning  of  observed  test  data.  As  a  result, 
one  may  well  reason  that  it  is  better  to  avoid  designing  for  operation  in  the 
boundary  lubrication  regime. 
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However,  devices  which  operate  in  the  boundary  lubricated  regime  are  many; 
the  lack  of  a  rigorous  mathematical  relationship  notwithstanding.  Obviously 
a  sound  equation  can  be  expected  to  improve  the  ability  to  successfully 
design  for  operation  under  boundary  lubricated  conditions.  For  this  effort 
the  ability  to  better  predict  scuffing  failure  load  of  relatively  simple 
sliding/rolling  contact  of  a  disc  pair  was  accepted  as  a  reasonable  objective. 
If  this  could  be  achieved,  perhaps  more  complex  devices  could  be  considered 
for  further  improvement  and  generalization  of  the  mode. 

2.  Development  of  the  Model 

Predicting  disc  load  at  scuffing  failure  has  already  been  addressed  using 
the  tools  of  thermodynamics,  chemistry,  physics,  etc.  in  various  combinations. 
The  approach  selected  for  model  development,  in  this  case,  was  based  simply 
upon  curve  fitting  of  experimental  data. 

Limited  resources  precluded  the  performance  of  new  disc  experiments 

specifically  designed  to  measure  and/or  control  all  the  suspected  parameters 

(4) 

of  importance.  Fortunately  the  appendix  of  an  Air  Force  reportv  '  contains 
some  excellent  experimental  data  including  many  of  the  parameters  of  interest. 
Unfortunately  some  of  the  parameters  reported  are  calculated  from  empirical 
relationships  or  relationships  having  unconfirmed  accuracy.  However,  these 
were  considered  for  inclusion  in  the  model  as  measured  values  were  not  readily 
available . 

The  cited  data  source  reports  on  280  disc  experiments  grouped  according 
to  the  type  of  lubricant  used,  method  of  lubricant  distribution,  etc.  Three 
of  these  basic  groups  were  selected  for  analysis.  As  a  result,  the  scope  was 
reduced  to  a  more  manageable  size  at  the  cost  of  generality. 

The  number  of  variables  considered  independent,  yet  related,  precluded 
utilizing  simple  methods  of  curve  fitting.  Fortunately  a  statistical  method 
known  as  multiple  regression,  a  type  of  response  surface  analysis,  has  been 
developed^ . 

To  use  this  method,  the  load  capacity,  W,  was  first  expressed  as: 

W  =  F  (G,  Vt,  Vs ,  f,  T0,  Ts,  Tc,  h.  A,  +)  (1) 
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Where : 


G  =  Lubricant  flow  rate,  gal /min 

Vf-  =  Sum  Velocity,  in/sec 

Vs  =  Sliding  Velocity,  in/sec 

f  =  Coefficient  of  Friction,  nondimensional 

T0  =  Conjunction  (Load  Zone)  Inlet  Oil  Temperature,  °R 

Ts  =  Average  Disc  Surface  Temperature,  °R 

Tc  =  Conjunction  Temperature,  °R 

h  =  Oil  Film  Thickness,  m  inches 

A  =  Specific  Film  Thickness.  The  ratio  of  lubricant  film 
thickness  to  the  composite  surface  roughness  of  the 
two  mating  surfaces,  nondimensional. 

<j>  =  Frictional  Power  Loss,  Btu/sec 

This  new  model  was  then  anticipated  to  assume  the  form: 

W  =  (GBl  Vt32  Vs33  f34  T03s  Ts3s  Tc3?  h38  A39  <j,3l0)C  (2) 

In  the  above,  the  p's  and  C  are  constants  to  be  determined. 

For  the  model  to  be  more  amenable  to  further  development  by  the  selected 
response  surface  analysis,  it  was  transformed  to: 

LnW  =  BiLnG  +  e2LnVt  +  B3LnVs  +  34Lnf  +  35Ln  (T0'+  459.7) 

+  36Ln  (Ts‘  +  459.7)  +  g7Ln  (Tc‘  +  459.7)  +  e8Lnh 

+  BgLnA  +  3ioLn<))  +  Bo  (3) 

Where 

Bo  =  Ln  C 
and 

T0‘,  Ts',  and  Tc'  can  now  be  inserted  in  °F. 
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Even  with  the  reduced  scope,  the  task  of  curve  fitting  90  data  sets  into 
the  above  is  overwhelming.  Fortunately  there  are  computer  library  routines 
which  allow  the  task  to  be  accomplished  with  reasonable  ease.  One  of  these 
routines  identified  as  REGRESS  is  available  in  University  of  Dayton's 
computer. 

Equation  3,  along  with  90  sets  of  data,  was  entered  into  the  computer  to 
solve  for  each  3  coefficient.  Values  for  all  coefficients  except  g9  were 
found.  Without  g9  the  F  ratio  statistic,  also  determined  by  REGRESS,  was 
found  to  have  the  value  of  843.284.  At  this  value  the  computer  associated 
alpha  was  zero  implying  the  probability  bound  had  been  reached.*  The  b's 
determined  are: 


Bo  =  8.41732 
3i  =  -0.532615 
=  -1.73055 
B3  =  0.052397 
B4  =  -0.0631378 
B5  =  2.97684 
B6  =  8.66441 
B7  =  -9.997 
B8 =  0.933837 
Bg  =  0.00000 
Bio  =  2 .07596 

Substituting  these  b's  into  Equation  3  and  then  taking  the  analog  results 


in: 


W  =  4524.76 


y  0.0523297  (y  '  +  459.7)2.97684  (y  1  +  459.7)8.66441 

q0. 0532615  y  1.73055  f 0 . 0631378  (y  '  +  459.7)9-997 
t  C 

^0. 933837  ^2.07596 


(4) 


*This  conclusion  is  due  to  the  limits  of  the  F  statistic  table  stored  in 
the  computer.  In  a  physical  sense  it  identifies  A  as  having  no  effect  on  W 
at  scoring  failure.  It  confirms  full  film  lubrication  does  not  exist  when 
scoring  occurs. 


8 


The  real  worth  of  Eq  4  was  assessed  by  investigating  the  proportion  of 
variability  in  W  that  it  could  be  explained.  This  was  accomplished  by  calcu¬ 
lating  the  quotient  of  the  regression  sum  of  squares  divided  by  the  total 
sum  of  squares.  For  Eq  4,  this  quotient,  the  correlation  coefficient,  was 
found  to  be  0.994771.  That  is,  98.9569  percent  of  the  variation  in  load, 

W,  can  be  explained  by  Equation  4. 

It  should  be  noted  this  statement  is  valid  for  mean  values  of  the 
dependent  as  well  as  the  independent  variables.  That  is,  Equation  4  cannot 
exactly  predict  W  as  all  the  variables  have  some  distribution  (assumed  normal) 
about  their  mean  value. 

3.  Comparison  of  the  Model  to  Experimental  Data 

Test  numbers  were  randomly  selected  from  the  90  data  sets  used  to  develop 
Equation  4.  The  independent  variables  from  15  of  these  sets  were  used  in 
Equation  4  to  compute  failure  load.  These  calculated  failure  loads  are  com¬ 
pared  to  experimentally  determined  values  in  Table  1.  The  error  values 
tabulated  are  based  upon: 

•  Assuming  the  measured  load  to  be  accurate  to  the  nearest  50  lb 
increment,  and 

•Adjusting  the  computed  load  to  the  nearest  50  lb  increment. 

The  mean  error  of  these  selected  data  sets  is  0.9%  with  a  standard  deviation 
of  8.0%.  Based  upon  the  assumption  that  the  variables  have  normal  distributions 
then,  by  the  use  of  t  statistics,  mean  error  is  predicted  to  lie  in  the  range 
from  -5.3%  to  7.0%  with  99%  confidence. 

It  was  recognized  that  the  worth  of  Eq  4  could  better  be  judged  on  its 
ability  to  predict  failure  loads  for  the  full  spectrum  of  data  reported  in 
Reference  4.  The  full  data  spectrum  (280  data  sets)  encompasses  a  considerable 
range  expansion  of  some  parameters  and  other  potentially  important  con- 
side  rati  ons . 

The  range  of  parameters  and  conditions  associated  with  the  basic  90  data 
sets  and  the  full  spectrum  of  groups  are  presented  in  Table  2  for  ease  of 
comparison. 
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TABLE  1 


COMPARISON  OF  COMPUTED  VERSUS  EXPERIMENTALLY 
DETERMINED  FAILURE  LOAD,  BASIC  GROUPS^ 

Computer  Load  Calculated  Load  Experimental  Load  Error 
Test  Nr  (b)  (Nearest  Pound)  (Nearest  50  Pounds)  (Nearest  50  Pounds)  (Percent) 


54 

1010 

1000 

950 

+5.3 

57 

513 

500 

500 

0 

29 

1068 

1050 

1100 

-4.5 

92 

541 

550 

500 

10.0 

97 

553 

550 

600 

-8.3 

84 

857 

850 

800 

6.2 

43 

2240 

2250 

2250 

0 

63 

648 

650 

650 

0 

71 

2992 

3000 

3300 

-9.1 

31 

1162 

1150 

1100 

4.5 

96 

653 

650 

750 

13.3 

20 

3362 

3350 

3450 

-2.9 

72 

1436 

1450 

1550 

-6.4 

94 

470 

450 

500 

-10.0 

99 

1155 

1150 

1000 

15.0 

(a) 

The  90  data  sets  form  the 

Basic  Groups. 

These  sets  were 

used  in  the 

development  of  Equation  4. 

(b) 

These  Test  Nrs  correspond 

to  those  used 

in  Reference  4, 

Appendix  D. 
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TABLE  2 


COMPARISON  OF  PARAMETER  RANGE  AND 
CONDITIONS,  BASIC  GROUPS  VERSUS  FULL  SPECTRUM 


Parameter 

or  Condition  Basic  Groups  All  Groups 


Disc  Material 

AMS  6265 

AMS  6263,  AMS  6475 

Lubricant 

MIL-L-7808G 

MIL-L-7808G  &  H 

Disc  Interacting 
Frequencies 

High  Speed  Disc  24  to 

30  revolutions/contact; 

Low  Speed  Disc  9  to  13 
revol utions/contact. 

Same  as  Basic  Group. 

Also  High  Speed  Disc  2 
to  7  revolutions/contact; 
Low  Speed  Disc  1  to  3 
revol uti ons/contact. 

Disc  Surface  Composite 
Roughness 

6.0  ±  1 .4  y  in  AA 

3.5  ±  0.7  y  in  AA  and 

6.0  ±  1 .4  y  in  AA 

Oil  Supply  Method 

Conjunction  Inlet 
or  Conjunction  Exit 
w/Disc  Side 

Conjunction  Inlet, 
Conjunction  Exit,  or 
Conjunction  Exit  w/Disc 
Side 

Oil  Supply  Rate 

0.033  to  2.5  gal/min 
total 

0.033  to  100  gal /mi n 
total 

Oil  Supply  Jet  Temp 

1 90°F 

190  and  275°F 

Ambient  Pressure 
(Disc  Environment) 

14.7 

14.7  psia  and  22.2  psia 

Sum  Velocity 

600  and  1200  inch/sec 

600,1200,2000,  and 

3000  inch/sec 

Sliding  Velocity 

178,  237,  292,  356,  474, 
and  583  inch/sec 

178,  200,  237,  240,  292, 
300,  356,  474,  583  inch/ 
sec 
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Test  sets  rrom  the  full  data  spectrum  were  randomly  chosen.  The 
independent  parameters  were  used  in  Equation  4  to  estimate  the  load  at 
failure.  These  computed  values  are  compared  to  the  experimentally 
determined  values  in  Table  3.  The  error  values  tabulated  are  based  upon: 

•  Assuming  the  Measured  Load  to  be  accurate  to  the  nearest 
50  lb  increment,  and 

•  Adjusting  the  Computed  Load  to  the  nearest  50  lb  increment. 

The  mean  error  of  this  sample  of  the  total  population  is  -2.6%  with  a 
standard  deviation  of  7.9%.  Based  upon  the  assumption  that  the  variables 
have  normal  distributions  then,  by  the  use  of  t  statistics,  the  mean  error 
is  predicted  to  lie  in  the  range  from  -8.6%  to  3.5%  with  99%  confidence. 
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TABLE.  3 


COMPARISON  OF  COMPUTER  VERSUS  EXPERIMENTALLY 
DETERMINED  FAILURE  LOAD,  ALL  GROUPS ^ 

Computed  Load  Computed  Load  Experimental  Load  Error 

Test  Nr  (b)  (Nearest  Pound)  (Nearest  50  Pounds)  (Nearest  50  Pounds)  (Percent) 


150 

773 

750 

750 

0 

2 

2250 

2250 

2300 

-2.2 

241 

1359 

1350 

1300 

3.8 

182 

1000 

1000 

950 

5.3 

239 

520 

500 

600 

-16.7 

213 

2677 

2700 

2900 

-6.9 

94 

470 

450 

500 

-10.0 

73 

3137 

3150 

3400 

-7.4 

180 

966 

950 

1000 

-5.0 

7U 

2879 

2900 

3150 

-7.9 

99 

1155 

1150 

1000 

15.0 

200 

2947 

2950 

3300 

-10.6 

186 

1791 

1800 

1800 

0 

43 

2240 

2250 

2250 

0 

137 

2336 

2350 

2250 

4.4 

a.  Two  hundred  eighty  data  sets;  all  the  available  data  contained  in 
Ref  4,  Appendix  D. 

b.  These  test  numbers  correspond  to  those  used  in  Ref  4,  Appendix  D- 
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SECTION  III 
CONCLUSIONS 


Equation  4  provides  a  reasonably  accurate  means  for  the  determination 
of  scuffing  limited  failure  load  on  the  AFAPL  Disc  Tester. 

Inspection  of  Equation  4  predicts  the  parameters  of  greatest  influence 
on  load  capacity,  as  determined  by  the  AFAPL  Disc  Tester  to  be: 

Tc  =  The  conjunction  temperature 

Ts  =  The  average  disc  surface  temperature 

The  next  greatest  influence  is  expected  from: 

T0  =  The  load  zone  temperature 

=  The  frictional  power  loss 

Vj.  =  The  sum  velocity  of  the  surface 

The  technique  of  multiple  regression  holds  considerable  promise  as  a  means  of 
developing  a  more  general  relationship  describing  disc  and  perhaps  gear  scuffing 
failure  limited  loads. 

Multiple  regression  analysis  is  anticipated  to  be  a  useful  tool  for 
investigating  other  complex  phenomena.  It  should  be  capable  of: 

•Providing  interim  empirical  relationships; 

•  Suggesting  the  most  promising  approaches  to  develop  more  rigorous 
relationships  from  first  principles;  and 

•  Identifying  significant  system  parameters,  with  their  relative 
strength  of  influence.  (Monitoring  of  these  system  parameters  may  provide 
the  needed  information  to  support  "on  condition  maintenance"  concepts.) 
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SECTION  IV 


RECOMMENDATIONS 

A  Complete  multiple  regression  analysis  of  data  from  the  AFAPL  Disc  Tester 
is  warranted.  The  analysis  should  result  in  a  more  generalized  relationship  with 
improved  accuracy.  This  analysis  should  be  followed  by  a  new  experimental 
program  with  this  machine. 

Any  new  experimental  program  should  maximize  the  parameters  measured  and/or 
controlled.  The  use  of  calculated  parameter  values  such  as  the  coefficient  of 
friction  should  be  avoided  where  possible.  Again,  an  improvement  in  accuracy  is 
anti  ci  pa  ted. 

Multiple  regression  analysis  should  be  investigated  as  a  means  of  establishing 
relationships  from  readily  measured  engine  parameters.  These  relationships 
(monitoring  engine  health)  may  be  combined  electronically  to  annunciate  not  only 
a  need  for  maintenance  but  what  the  problem  may  be. 
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